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We present results of the first-principles calculation of Cs dipole static polarizabilities for the Ns 
{N = G - 12), Npj (iV = 6 - 10), and Ndj {N ^ 5 ~ 10) states using the relativistic all-order 
method. In our implementation of the all-order method, single and double excitations of Dirac-Fock 
wave functions are included to all orders in perturbation theory. Additional calculations are carried 
out for the dominant terms and the uncertainties of our final values are estimated for all states. A 
comprehensive review of the existing theoretical and experimental studies of the Cs polarizabilities 
is also carried out. Our results are compared with other values where they are available. These 
calculations provide a theoretical benchmark for a large number of Cs polarizabilities. 



I. INTRODUCTION 

Atomic and molecular polarizabilities are of great interest because various properties of matter 
can be expressed in terms of multipole moments and polarizabilities of the atoms or molecules in 
the system. Polarizabilities describe the response of a system to external electric fields. Therefore, 
atomic polarizabilities reflect the atomic structure and can be used to probe correlation and relativistic 
effects. On the other hand, atomic polarizabilities determine the long-range van der Waals interactions 
between the atoms and are used in describing atomic scattering processes. Atomic polarizabilities play 
an important role in high Rydberg state spectroscopy (e.g. [1|). The study of the alkali- metal atoms 
are of particular interest because they allow for very accurate comparison between the experiment and 
theory. The conclusions reached from such studies may also be useful for the understanding of more 
complicated systems. Cs is also of particular interest owing to the study of the parity nonconservation 
(PNC), designed to test the standard model of the electroweak interaction and to set limits on its 
possible extensions as well as to infer nuclear anapole moments. The highest accuracy of such an 
experiment was reached for cesium, where measurements of PNC amplitudes have reached an accuracy 
of 0.4% 0. To make meaningful tests of the standard model, high-precision calculations of the PNC 
amplitudes must be carried out at a similar level of accuracy, and the uncertainty of the calculation 
has to be established. 

In this work, we carry out a systematic study of a large number of Cs polarizabilites in order to 
provide recommended values for the Ns {N = 6 — 12), Npj {N = 6 — 10), and Ndj {N = 5 — 10) states 
and evaluate their uncertainties. The best-set values for the 91 electric-dipole matrix elements used in 
our calculations are also provided with their uncertainties. These data are also useful for a number of 
other applications. 



A. Experimental methods and studies of the atomic polarizabilities 

In this section, we provide a summary of a variety of methods used to measure the atomic polariz- 
abilities as well as describe the development in the experimental measurements of the electric-dipole 
polarizability of cesium ground state. 
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In 2005, Gould and Miller 0] wrote a comprehensive review of the experimental methods to de- 
termine the static electric-dipole polarizabilities. Miller and Bederson's earlier review from 1988 [1] 
concentrated on the bulk polarizability measurements and the atomic beam methods. Average bulk 
ground state static polarizabilities are measured by determining the dielectric constant of an atomic 
or molecular gas. The bulk dynamic polarizabilities are determined by measuring the refractive index 
of the gas, see The bulk methods are very accurate, but their limitation lies in the need to deal 
with atoms or molecules that are stable and gaseous at room temperature and the fact that the effect 
of the excited states can not be accounted for. 

In 1974, Molof et al. used the E-H-gradient balance technique to measure the static electric- 
dipole polarizabilities of alkali-metal atoms. They obtained the value (59.6 ± 1.2) 10^^^ cm'^ for 
electric-dipole polarizability of the ground state of cesium. Hall and Zorn measured the value (63.3 
± 4.6) 10"^"* cm'^ for the electric-dipole polarizability of the ground state of cesium. They used the 
deflection of a velocity-selected atomic beam in inhomogeneous electric field. The technique is based 
on the fact that the deflection experienced by atoms moving through a region with known transverse 
electric field gradient is proportional to the dipole polarizability of the atoms. An important detail of 
this technique is that the precision with which the velocity of the atoms is known puts a limitation on 
the precision of the experiment. The short interaction time in the case of high velocity which leads to 
small deflection of the beam places another limitation on the accuracy of this method. 

In 1995, Ekstrom et al. f?! designed an atomic interference experiment that allowed them to measure 
the ground state energy shifts with spectroscopic precision and determine the ground state dipole po- 
larizability. In 2003, Amini and Gould designed an experiment that avoids the problems associated 
with the measuring the deflection of a thermal beam in transverse electric-field gradient. They measure 
the effect of the electric-field gradient on the longitudinal velocity of the beam of cesium atoms in a 
magneto-optical trap (MOT). The cesium 6s scalar dipole polarizability is found from the time-of-flight 
of laser cooled and launched cesium atoms traveling through an electric field. The cited value is (59.42 
± 0.08) 10^^** cm'^. This is the most precise measurement of the ground state polarizability at this 
time. 

Another group of experiments allows to infer the atomic polarizabilities by measuring the Stark 
shift of the cesium lines, e.g. 0. In 1988, Tanner and Wieman [l3| measured the Stark shift in 
the 6si/2 — * 6P3/2 transition in Cs. The dc Stark shift of the cesium Dl line has been has been 
measured to 0.01% in Ref. [ll[. The authors of this work [ll| noted that it was the most precise Stark 
shift measurement ever reported. The Stark shifts of the 6^3/2 — (10 — 13)s states in cesium were 
measured in Ref. [l2|. The Stark shifts of cesium IID states were measured with high precision by 
van Wijngaarden and Li in 1997 fl^ using an electro-optically modulated laser beam. The authors 
note that the tensor polarizabilities reported in that work [13] were the most accurate yet determined 
for any atomic state. The dc Stark shift of the 6s — 7s transition in atomic cesium was measured with 



14 1 using laser spectroscopy. The result of this experiment disagrees with a 
15l | but was within 0.3% of the value predicted by the ah initio calculations 



high precision in 1999 
previous measurement 

jl6l . [l7j removing the largest at that time outstanding disagreement between experiment and ah initio 
theory of low-lying states in atomic cesium. 

The atomic polarizabilities can be derived from measurements of the lifetimes of the corresponding 
levels. The contribution of the core electrons to the polarizability of the alkali atoms is small. Then, 
the main contribution to the ground s-state polarizability comes from the first lo w-ly ing excited P- 
states, i.e. dominant electric-dipole matrix elements are (A^s|D|npi/2,3/2); see [3 [13 ^ detailed 
discussion and comparison of lifetime and polarizability measurements in cesium. 

A large group of experiments makes use of the level-crossing of some hyperfine atomic levels at 
finite electric field. The first observation of the purely electric field level-crossing was reported in 1966 
[20j . This type of measurements allows for experimental determination of the excited states tensor 
polarizabilities. Recent cesium measurements were reported by Auzinsh et al. [21I [2^. 
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B. Theoretical studies of cesium polarizabilities 

Since the alkali-metal atoms are monovalent systems, they represent an excellent opportunity to 
study the correlation effects. Heavy alkalis are of particular interest owing to the application to the 
study of fundamental symmetries. The polarizability of the alkali-metal atoms are essentially the same 
as the valence polarizability as the contribution of the ionic core was determined to be small [1^, HI] • 
We summarize the theoretical studies of Cs polarizabilities below. 

In his seminal paper [25| . Dalgarno summarized the mathematical foundation of the theory of the 
atomic perturbation and discussed the methods of calculating the atomic polarizabilities and shielding 
factors. The polarizability of the cesium atom obtained by Dalgarno and Kingston [2^ using the 
oscillator-strength formula was (53.7 ± 5.4) 10"^'* cm^. 

According to the oscillator-strength formula, the knowledge of the (reduced) electric-dipole matrix 
elements is crucial for calculation of the atomic polarizabilities. The reduced matrix elements can be 
computed in a number of approximations. Variety of theoretical methods are used, such as third-order 
many-body perturbation theory, multiconfiguration Hartree-Fock (MCHF), configuration interaction 
(CI) method, coupled-cluster (CC) method, and it relativistic linearized version referred to as the 
all-order method as well as others. 

In 1970, Sternheimer (27| used the Hartree-Fock wave functions to compute the quadrupole polar- 
izability of some ions and alkali atoms. The cesium ground state value was calculated to be 71.31 
10^^"* cm'^. In 1971, Schmieder et al. [11] calculated the scalar and quadrupole polarizabilities of 
cesium states in the second order perturbation theory. The work by Kello et al. [l^ contains 
a detailed investigation of the contracted Gaussian basis sets in the calculation of the electric-dipole 
polarizabilities of alkali-metal atoms. The calculations are performed using the complete-active-space 
self-consistent field and second order perturbation theory, CASSCF and CASPT2. Another group of 
Gaussian basis set methods use relativistic pseudopotentials (see (s^j and the references there). Pseu- 
dopotential methods replace the core electrons by an effective, pseudopotential. The core polarization 
and the relativistic effects are incorporated as well. The Douglas-Kroll relativistic CCSD(T) method 
with the optimal basis set gives 58.09 10~^* cm'^ for the cesium ground state dipole polarizability. 

Extensive calculation of the polarizabilities of cesium Ns, Np, Nd, and nF states was carried out 
by van Wijngaarden and Li [31| using the Coulomb approximation. They also provided extensive 
comparison with other theoretical and experimental values. 

Patil and Tang [33| computed the multipolar polarizabilities, aq, with q = 1,2,..., 12, for the 
alkali isoelectronic sequences. The ground state wave functions were taken to be the asymptotically 
correct wave functions, i.e. the two leading terms in the asymptotic expansion of the wave function 
are retained. The excited states are taken to be the Coulomb wave functions with a correction that 
makes sure the experimental energies of the low-lying states are reproduced correctly. The ground 
state electric-dipole polarizability of cesium was found to be 60.6 10^^^ cm"^. 

The relativistic linearized version of the coupled-cluster singles-doubles method, i.e. all-order SD 
method, was used in [2^ [23| to calculate the static dipole polarizabilities of the alkali-metal atoms. 
This method is discussed in more details in Section jllll The value obtained for the ground state static 
dipole polarizability is 59.3(3) 10"^** cm'^ [2^. In [S^, Porsev and Derevianko computed the ground 
state quadrupole and octupole polarizabilities of the alkali-metal atoms using the relativistic MBPT. 

In 2004, Safronova and Clark [l^ pointed out the inconsistencies between the lifetime and polar- 
izability measurements in cesium. The dominant contribution to the 6p scalar polarizability comes 
from the 5c? — 6p matrix elements. This allows for a check of the accuracy of the matrix elements. 
The authors compare the values of the 6p polarizabilities obtained using the relativistic all-order SD 
method and using the values of the matrix elements derived from the 5d lifetime experiment [3^ . 
They point out that the theoretical all-order values yield a polarizability result in agreement with the 
polarizability measurements [13, [ii| but not with the lifetime measurements [IJl . 

In a recent work, Gunawardena et al. (ssj report a precise measurement of the dc Stark shift of 
the 6s — > 8s transition in atomic cesium. The experiment makes use of the Doppler-free two-photon 
absorption measurement. The value of the static polarizability of 8s state in cesium, extracted from 
the experiment, is 38 060±250 Og. The authors present a theoretical value of 38 260±290 Gq. The 
theoretical value of the polarizability is calculated using the relativistic all-order SD method. 
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II. METHOD 

The energy shift of the \jjm) atomic level in a dc electric field E = £z is given by 



AE - --a^.rnS^, (1) 



where a^jm defines the static polarizability of the corresponding atomic state |7jm). The scalar and 
tensor static polarizabilities ao^-yjmj and a2,'yjm are defined as 

3m^-j(j + l) 

We omit index jjm in the text below and refer to scalar and tensor static polarizabilities as ao and 
a2, respectively. 

We separate the calculation of the scalar static polarizability into the calculation of the polarizability 
of the ionic core and the valence polarizability. The random-phasc-approximation (RPA) calculation 
of the Cs core polarizability was carried out in Ref. [36] and yielded the value IS.Soq, where cq is the 
Bohr radius. Based on the evaluation of the accuracy of RPA approximation for the polarizabilities 
of the noble gases, this value is accurate to at least 5%. The core polarizability is small even in 
comparison with the valence polarizabilities of the lowest states such as 6p and 5d. It is negligible 
for the higher states. For example, core polarizability contributes only 4% to the total value of the 
ground state polarizability and only 1% to the 6pi/2 polarizability. Therefore, the RPA value of the 
core polarizability is sufficiently accurate for this work. The separation of the scalar polarizability to 
the core and valence parts also produces a compensation term that accounts for the Pauli exclusion 
principle, i.e. in Cs it subtracts 1/2 of the core polarizability contribution associated with the excitation 
to the valence shell. This term is only 2% of the core contribution [23] even for the ground state and, 
therefore, below the estimated uncertainty of the core term itself. It is negligible for all other states. 

The valence scalar and tensor static polarizabilities of the atomic state Ijj) are expressed in lowest 
order as sums over unperturbed intermediate states {(ijp) of parity opposite to that of the state \ 

ao = ^ 



Oi2 



5j(2j - 1) r j 1 3p \ \{ij\\Dmp)\ 



6(2j + 3)(2j + l)(j + l)^^ > 11 J 2/ Ep~E, ' 



where (7j11£'11/3j^) is the reduced electric-dipole matrix element defined as 

{^MD.Wj'm') = i-iy'-^' J I ) {uWDWif), (5) 

and the Dq is the corresponding component of the electric-dipole operator in spherical coordinates. 

The sums over states (3 in Eqs. (|3l4p separate into the two or three sums over the principal quan- 
tum number for each type of the allowed electric-dipole transitions for Cs calculation. The al- 
lowed values of (3 are the following: [3 = 7 '^■^3/2 for states, (3 = ns,nd^^2 for npi/2 
states, l3 = ns,nd^/2,nd^/2 for the np^/2 states, (3 = npi/2, ?^P3/2j "■/5/2 for the nd'j,/2 states, and 
(3 = np3/2, n/5/2, n/7/2 for the 71^5/2 states, here, n and A'^ are principal quantum numbers. We note 
that there is no tensor contribution to the polarizability of the Ns and Npi/2 states. 

In order to evaluate the sums over the principle quantum numbers n, we carry out all calculations 
in a finite B-spline basis set ^3] constrained to a large spherical cavity and defined on a non-linear 
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grid. Therefore, the sums in Eqs. (|3l4p range over the basis set states. In order to study such highly- 
excited states as 12s and 10c?, we needed to use a very large cavity, R — 220 a.u, to ensure that the 
corresponding wave functions fit inside the cavity. As a result, we had to use a large number of basis 
set functions, Nb = 70, in order to correctly represent the properties of both highly-excited states and 
the lowest states. We verified that the basis set energies match the Dirac-Fock energies for all of the 
states considered in this work. The order of splines was taken to be fc = 8. More calculation details 
associated with the use of the finite basis set are described in Section IIIII 

The sums over the principal quantum numbers in Eqs. pi4p converge very rapidly, with the exception 
of the sums involving the 5c?3/2 — «/5/2 a-nd 5^5/2 — /7/2 transitions which we will discuss separately 
in Sectior flVl In fact, each of the sum over the principal quantum number is dominated by one or 
two terms that correspond to the lowest possible values of the denominator — Ej in Eqs. pi4p . A 
small number of other terms may also be significant for the precise calculation for some states. As a 
result, only a few terms from each sum have to be calculated accurately, allowing us to separate the 
calculation of the valence scalar and tensor polarizabilities into the two parts, the main term containing 
all significant contributions and the tail: 

„ „ main i ^ tail /a\ 

"0,2 ="0,2 +"0,2- (6) 

The separation of a into the main and tail parts is done independently for each of the two or three 
sums over the principal quantum number n contributing to the polarizability of the specific state: 

"main Nb 
n nn ?lmain + l 

where the uq is the lowest possible value of principal quantum number for the particular set of (3 
states, rimain is the cut-off that we select for the separation of the main and tail terms, and Nb is 
the number of the finite basis set orbitals set to 70 in the present work. In general, rimain may be 
selected differently for the specific state 7 and each (3 sum, but we chose to use the same nmain for 
all of the states considered in this work. We use rtmain = 12 for (3 = ns, «Pi/2, "■P3/2J "•main — 10 for 
P = ndi/2,nd2,/2, and rimain = 8 for j3 = ?^/5/2, 't-c?7/2j respectively. The only exception is the addition 
of the 9/ contributions to the calculation of the lOd polarizabilites. 

Such high value of the cut-off principal quantum number also reduced the tail contribution and 
improved the accuracy of our calculations. The remaining tail contributions are evaluated using in the 
Dirac-Fock (DF) approximation, i.e. both energies and El matrix elements were calculated in the DF 
approximation. 

In summary, we reduce the calculation of the polarizabilities to the calculation of the electric-dipole 
reduced matrix elements required for the evaluation of the main terms for each state. We use the 
experimental energies from Refs. [H, [s^ H^l in the calculation of the main terms. Owing to the 
large number of states considered in this work, 317 transitions contribute to the main term and 102 
transitions give dominant contributions. We calculated all of the 317 electric-dipole matrix elements 
using the relativistic all-order method and conducted additional calculation for the 102 transitions 
that involved the evaluation of the largest missing corrections and evaluation of the uncertainty of the 
final values. The calculation of the matrix elements is described in the next section. 

III. CALCULATION OF THE El MATRIX ELEMENTS 

We carry out the calculation of the electric-dipole reduced matrix elements using the relativistic SD 
all-order method where single and double excitations of the (frozen-core) Dirac-Fock wave function are 
included to all orders in perturbation theory [l^, EH ■ Triple excitations are also partially included for 
selected cases. The relativistic SD all-order method is a linearized coupled-cluster method restricted 
to single and double excitations. A comprehensive review of the coupled-cluster method and its 
applications in quantum chemistry is given in Ref. [4^ . In the coupled-cluster method, the exact wave 
function of the monovalent atom in a state v is represented as 



1^-,) =exp(5)|$,), 



(8) 
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where is the lowest-order atomic wave function for the state v, which taken to be a frozen core 
Dirac-Fock (DF) wave function in our calculations. The cluster operator S is expressed as a sum of 
n-particle excitations Sn of the lowest-order wave function 

S = Si + S2 + S3 + ---. (9) 

The exponential function in Eq. ([8]) can be expanded to give 

|*.>-(i + ^ + ^^' + ---)l'i>.)- (10) 

In the linearized single-double (SD) coupled-cluster method, only terms that are linear in the Si 
remain and all other terms, for example Si x 52 are omitted , i.e. 

= (1 + 51-1-52 + ■••) I*.). (11) 

The contributions from the non- linear terms were recently investigated in Refs. [43', 'ii', '45*1 . We refer 
the reader to Ref. 45] for a complete list of the non-linear terms and detailed investigation of their 
contributions to the alkali-metal atom properties. The computational complexity of the calculations 
increases dramatically beyond the double excitations term 52 , and we include triple excitations partially 
in some of the calculations using a perturbative approach. We note that in this work very large 
{Nb — 70 for each partial wave) complete basis set is used to correctly reproduce necessary atomic 
properties for highly excited atomic states, requiring significant computational resources for the SD 
all-order calculations. 

The expression for the single excitations is given by 

5l = Pmaalntta + ^ Prnval^ay, (12) 
ma m^v 

where the first term corresponds to single core excitations and the second term corresponds to single 
valence excitations. The expansion coefficients pma and pmv are referred to as single core and valence 
excitation coefficients, and the a| and are creation and annihilation operators for an electron in a 
state i. We use the letters from the beginning of the alphabet a,b,... to designate core orbitals and 
letters from the middle of the alphabet, m, n, ... to designate excited states. For Cs, we include all 17 
a = Is, 5^3/ 2 core shells in our calculations. 
The double excitation term is given by 

S2 — ^ Pmnabaln'A.O'baa + PmnvbO'mO'n^bav i (13) 

mnab mnb 

and the quantities Pmnab and Pmnva are referred to as double core and valence excitation coefficients, 
respectively. 

Therefore, the atomic wave function in the all-order SD method [4l| is expressed via the single and 
double excitation coefficients as follows: 

1^'^°) = ( 1 + X! P"^o.alnaa + Pmriabai^alataa 

\ ma mnab 

+ ^ PrnvOl^ay + ^ PmnvaalnO^aaa^ \'^v)- (14) 
m^v mna j 

The equations for the excitations coefficients and the correlation energy are derived by substituting 
the SD all-order wave function given by the Eq. (|14p into the Schrodinger equation 



(15) 
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where the Hamihonian H is the relativistic no-pair Hamihonian [46,] , which can be written in second- 
quantized form as 



H '^'"1'^* 9i]kia\a]aiak - ^ Uija\aj, (16) 

i ijkl ij 

where are the one-body DF energies for tlie state i, Uij is taken to be frozen-core DF potential in 
our calcuatHon, and gyfe; are the two-body Coulomb integrals: 

g,^,i= J d'r J d'r'^plir)^]ir')y^^Mr)Mr')- (17) 
For example, the equation for the single valence excitation coefficients pmv is given by 

i^v ^~ ^ ^v) Prav ^ ^ 9rahvnPnb ^ ^ 9mbnr Pnrvb ^ ^ QhcvnPranhcj (1^) 

hn bnr ben 

where SEy is the correlation correction to the valence energy for the state v given in terms of the 
excitation coefficients by 

— ^ ^ 9vavmPma ^ ^ 9abv7nPnivab "t" ^ ^ 9vbm7i Pjnnvb • (1"^) 
ma niab mna 

We use the designation pmnab = Pmnab-Pmnba and gmnab = 9mnab-9mnba in the equations above. The 
sum over the magnetic quantum numbers is carried out analytically and the resulting equations are 
solved iteratively for the excitation coefficients until the correlation energy converges. The excitation 
coefficients are then used for the calculation of tlic matrix elements as described below. 

In general, the one-body operator Z can be written in second quantization as Z = Zija\aj. 
The expression for SD matrix elements of operator Z is obtained by substituting the SD wave function 
given by Eq. (fT4| into the expression 



^{^■^\^y,){^y\^y) 

The resulting SD matrix element is given by 

^ ^ Zy,y + + . . . + 



(21) 



where Zyjy is the DF matrix element, terms Z*^°\ . . . Z^*^ are linear or quadratic functions of the 
excitation coefficients, and Ny and iV^j are normalization terms that are quadratic functions of the 
excitation coefficients. For most of the dominant transitions in our polarizability calculations, a single 
term 

Z^ ^ ^ ^ ZwraP^nv ~i~ ^ ^ Zj^y Pj^inj (22) 
m m 

gives the dominant contribution. Two other terms, 

Z^ ^ — ^ ^ ZamPwinva ~t- ^ ^ ZijidPyfniyd (23) 
ma ma 

and 

Z^ ^ ^ ^ ZmnPmwPnv (24) 
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may be dominant for selected important transitions. We note that both Z^'^-' and Z'^ terms contain 
only single valence excitations coefficients. The complete expression for the matrix elements is given 
in Ref. ^41i] . 

All sums over the excited state in the formulas above range over the basis set states. We truncated 
last five basis set orbitals for each partial wave since their contributions is negligible, i.e. 65/70 orbitals 
are included for each partial wave. All partial waves are included up to Imax = 6, and orbitals with 
2 = Z + 1/2 and j = / — 1/2 are considered separately since it is intrinsically relativistic calculation. The 
basis set is numerically stable, i.e. the increase of the number of the basis set orbitals does not change 
the results within the current accuracy. The numerical uncertainty associated with the truncation of 
the number of the partial waves at Imax = 6 is also negligible. We estimated the contribution from 
higher partial waves to be 0.1% for the 6s — 6pj transitions. The evaluation of the higher partial wave 
contribution is carried out by conducting the third-order perturbation theory calculation carried out 
as described in Ref. [l^l with the same basis set and with higher number of the partial waves. We 
also verified that the use of the very large cavity did not affect the numerical accuracy of the atomic 
properties of the lower states by comparing the present results for the transitions between the lower 
states with all-order calculation carried out with the small cavity appropriate for the lower states. We 
note that large basis set size, Nb — 70, is necessary to reproduce the atomic properties correctly with 
such a large cavity. We found that the accuracy of the Nb — 50 B-spline basis set is not sufficient for 
such calculation. 

As noted above, we have identified that the correlation correction for most of the dominant tran- 
sitions in our polarizability calculation is essentially determined by a single term, Z'^^^ that contains 
only single valence excitations. This term mostly corresponds to the Brueckner orbital correction as 
classified in the Ref. [l^]. It is established ([l6l. [23 , l48l . l49j that it can be corrected by scaling the 
single excitation coefficients using the ratio of the "experimental" and theoretical correlation energies 
and redoing the matrix element calculation with modified excitation coefficients. The "experimental" 
correlation energies are determined as the differences of the experimental data and our lowest-order 
DF values. We carry out such scaling where appropriate and determine the uncertainty of our calcu- 
lation of the matrix elements as the difference between the ah initio and scaled data. In certain cases 
where this correction is particularly large, we also carried out ah initio calculation of the limited triple 
excitations and conducted more accurate study of the uncertainty of the calculations. The limited in- 
clusion of the triples was also aimed at correcting the pmv excitation coefficients. Such calculations are 
described in detail in Refs. [13,111, SI] and references therein. We note that term Z'^'^^ is also corrected 
by scaling as it contains only valence single excitation coefficients. We note that the scaling procedure 
allows to place an uncertainty on our theoretical data that is not derived from the comparison with 
the experiment. Our results are summarized in the next section. 



IV. RESULTS AND DISCUSSION 



The results of the all-order calculation of the matrix elements are summarized in Table [D Owing to 
the very large number of the transitions involved in this calculation, we only listed the transitions that 
give dominant contributions to the polarizabilities of the states considered in this work. In order to 
provide a best set of known data for these transitions, we replaced all-order theoretical values by the 
experimental ones where high-precision values are available. The experimental values from Ref. [50| 
are used for the 6s — 6pj transitions, the values for the 6pj — 7s transitions are derived from the 7s 
lifetime measurement in [5lj , and the 6s — 7pj values are experimental values from [5l| . The 7s — 7pj 
values were derived from the 7s — 6s Stark shift measurement [T^ . We are quoting these 7s — 7pj values 
in the present table as the most accurate values available, and we used them in the calculation of the 
polarizabilities to provide recommended values for these states. However, we used our theoretical 
values in the calculation of the 7s polarizabilites for the evaluation of the accuracy of our calculation. 
Otherwise, the comparison of the 7s values with the experiment would have provided no information 
as we would have expected near exact agreement. Our theoretical values, 10.31(4) and 14.32(6), are 
in perfect agreement with values derived from the Stark shifts. 

The values for the 5(i — np and 6(i — 6p transitions are taken from the study of the inconsistencies 
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TABLE I: Absolute values of the selected reduced electric-dipole matrix elements El in Cs and estimates 
of their uncertainties. Unless otherwise noted, these are all-order SD scaled values, including values from 
Refs. [2^. |35|. "Experimental values, Ref. [SO^, 'SD all-order scaled values, previously published in Ref. [s^, 
"^experimental values from Ref. 51], ''derived from the 7s — 6s Stark shift value in Ref. (23), ''all-order values, 
Ref. [3, ■'^SD all-order scaled values, previously published in Ref. [22]. Units: eao. 



Transition 


El 


Transition 


El 


Transition 


El 


6s — 6jDi 10 


4.489(7)" 


8s — 7pi /o 


9.313(65)* 


10s - 9pi/2 


24.50(10) 


6s — 7pi 10 


0.276(2)'= 


8s — 8i)i /o 


17.78(7) ^ 


10s - lOpi /9 


38.31(10) 


6s — 6ps/9 


6.324(7)" 


8s — 7?33/9 


14.07(7) ^ 


lOs — 9P3 /9 


36.69(10) 


6s ~ 7p'^/9 


0.586(5)'= 


8s — 8p3/9 


24.56(10)'' 


10s - 10p'-,/9 


52.67(16) 


7s - 6pi/2 


4.236(21)'= 


9s - 8pi/2 


16.06(8) 


lis - lOpi/2 


34.64(12) 


7s — 7pi/2 




9s — 9pi/2 


z / .lU(^o ) 


1 1 CI 1 1 n 

lis — llpi/2 


CI A0^^ ^ \ 
oi.4zt^ i i ) 


7s — 6p3/2 


6.473(32)'= 


9s — 8p3/9 


24.12(8) 


lis — 10p3/9 


51.77(12) 


7s — 7p3/2 


14.320(14)'* 


9s — 9p3/9 


37.33(13) 


lis — llp3/9 


70.58(19) 


12s - llpi/2 


46.49(15) 


50(3/2 - 6pi/2 


7.06(16)" 


54/2 - 6p3/2 


9.66(20)'= 


12s - 12pi/2 


66.43(13) 


5^3/2 - 6p3/2 


3.19(8)'= 


54/2 - 4/5/2 


1.93(30) 


12s - llp3/2 


69.37(15) 


5^3/2 - 4/5/2 


7.1(5) 


54/2 - 4/7/2 


8.6(6) 


i2S — i2p3/2 


\)i.l[Z ) 










60(3/2 - 6pi/2 


4.15(20)" 


60(5/2 - 6p3/2 


6.01(26)'= 


74/2 - 7pi/2 


6.56(2)^ 


6^3/2 - 7pi/2 


18.0(2) 


64/2 - 7p3/2 


24.4(3) 


74/2 - 8pi/2 


32.0(2)^ 


60(3/2 " 6p3/2 


2.05(9)'' 


64/2 - 4/5/2 


6.60(5) 


74/2 - 7p3/2 


3.32(2)^' 


6^3/2 - 7p3/2 


8.07(11) 


60(5/2 - 5/5/2 


1.11(15) 


74/2 - 8p3/2 


14.35(8)^ 


60(3/2 - 4/5/2 


24.6(2) 


64/2 - 4/7/2 


29.5(2) 


74/2 — 4/5/2 


13.0(2)^ 


00,3/2 " 0/5/2 


o.9(D ) 


605/2 — 0/7/2 


4.9d(D / ) 


'4/2 — 0/5/2 




70(5/2 - 7p3/2 


9.64(4)-^ 


84/2 - 8pi/2 


9.18(5) 


84/2 - 8p3/2 


13.65(7) 


70(5/2 - 8p3/2 


43.2(2)^ 


84/2 - 9pi/2 


49.3(2) 


84/2 - 9p3/2 


66.6(2) 


70(5/2 - 5/5/2 


11.66(7)-'' 


84/2 - 8p3/2 


4.71(01) 


84/2 — 5/5/2 


6.85(4) 


70(5/2 - 4/7/2 


15.3(2)^' 


84/2 - 9p3/2 


22.13(7) 


84/2 — 6/5/2 


17.54(8) 


70(5/2 — 5/7/2 


52.2(3)-^ 


84/2 — 5/5/2 


26.1(2) 


84/2 — 5/7/2 


30.6(2) 






84/2 ~ 6/5/2 


65.2(4) 


84/2 — 6/7/2 


78.4(4) 


9^3/2 - 9pi/2 


12.2(2) 


94/2 - 9p3/2 


18.3(2) 


104/2 - lOpi/2 


15.6(2) 


90(3/2 ~ lOpi/2 


70.0(2) 


94/2 - 10p3/2 


94.5(2) 


104/2 - llpi/2 


94.1(2) 


90(3/2 - 9p3/2 


6.33(6) 


94/2 - 7/5/2 


24.36(9) 


104/2 - 10p3/2 


8.16(7) 


90(3/2 " 10p3/2 


31.45(8) 


94/2 - 6/7/2 


49.3(3) 


104/2 - llp3/2 


42.30(9) 


9d3/2 - 6/5/2 


42.0(4) 


94/2 - 7/7/2 


108.9(4) 


104/2 - 7/5/2 


61.0(2) 


9d3/2 - 7/5/2 


90.5(4) 






104/2 - 8/5/2 


119.4(4) 


100(5/2 - 10p3/2 


23.5(3) 


104/2 - 8/5/2 


32.2(1) 


104/2 - 8/7/2 


143.8(5) 


100(5/2 - llp3/2 


127.1(3) 


104/2 - 7/7/2 


71.7(3) 







in the lifetime and polarizability measurements in Cs [18| . The Sd — np values are scaled all-order 
values with the uncertainty defined as the spread of the ah initio values containing the partial triple 
excitations (SDpT) and scaled values based on both SD and SDpT calculations. The evaluation of the 
uncertainty of these values is described in [l^] . The 6p — Qd values are ah initio SDpT results with 
the uncertainty defined as the difference of the ah initio SD and SDpT calculations. We did not use 
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TABLE II: The contributions to the scalar polarizability for the 9s state in cesium. The corresponding energy 
differences and the absolute values of the lowest-order Z^^ and final all-order electric-dipole reduced matrix 
elements Z^^ are also listed. The energy differences are given in cm~^. Electric-dipole matrix elements are 
given in atomic units (eao), and polarizabilities are given in 10'^ af,, where ao is Bohr radius. 





B 


Z?^ 

/3,9s 


/3,9s 


J^8s 




main/ p \ 










n nn 




7Pl/2 


2.04 


1.96 


-5145 


-0.05 




8Pl/2 


-LU.OU 


1 fi Ofi 

-LU.UU 








9pi/2 


28.17 


27.10 


726 


74.0(4) 




lOpi/2 


2.67 


2.76 


1816 


0.31 




llPl/2 


1.01 


1.08 


2493 


0.03 




12pi/2 


0.56 


0.60 


2942 


0.01 


«*""(nPi/2) 










0.01 


a--"(nP3/2) 


6P3/2 


0.79 


0.77 


-15178 


0.00 




7p3/2 


2.86 


2.73 


-4964 


-0.11 




8p3/2 


24.31 


24.12 


-1119 


-38.0(3) 




9P3/2 


38.99 


37.33 


771 


132.3(9) 




10p3/2 


4.43 


4.61 


1843 


0.85 




llP3/2 


1.80 


1.93 


2510 


0.11 




12p3/2 


1.04 


1.13 


2954 


0.03 


a*""(nP3/2) 










0.04 


Total 










153.7(1.0) 



the measured lifetimes values [Mi owing to the inconsistencies of these values with the measured 
Stark shifts [Hi- The uncertainties of the 5^3/2 — 4/5/2 and 5^5/2 — 4/7/2 matrix elements are 
estimated as the differences of the SD scaled results and ah initio SDpT values that partially include 
triple excitations. 

The SD all-order values for the 8s — np and Idj — nlj transitions were previously published in 
Refs. [2^ d^, respectively. In summary, the uncertainties of our calculations are generally small, 
ranging from 0.2% to about 1%. The only exceptions are the transitions involving the 5d states and 
some of the transition from the 6d states. We refer the reader to Ref. [18i] for a detailed discussion 
of these transitions. We note that we may overestimate the uncertainty of our calculation for these 
transitions as our values for the Gpif2 and 6^3/2 polarizabilities are in excellent agreement with the 
experiment [lol.lTll|. As a result, the actual accuracy of our values of 5d polarizabilities may be actually 
higher than we estimated. 

As noted above, we used experimental energies for all of the main term calculations. Most of the 
energies values in this work are taken from the 1987 measurements by Weber and Sansonetti [39| and 
other values quoted in the same reference. The ionization potential value, required for the scaling 
procedure, is taken from the same work. The values of the several lower levels are taken from the 
NIST Handbook of Basic Atomic Spectroscopic Data [11]. The data for the np3/2 levels are obtained by 
combining the npi/2 values from [so'l and fine-structure intervals from [40| . The data for the remaining 
few levels not given in either |38i, [39] were taken from Ref. ■ Since the energy denominators in the 
polarizability calculation are small for some of the higher states (below 100 cm~^), we compiled the 
list of the most accurate known energies. As a result, the polarizability values quoted in this work 
for the 7d, 9d, and lOd states are slightly different from the ones quoted in Ref. 22] while the same 
matrix elements were used. We note that these differences are well within the uncertainties of the 
polarizability values. The uncertainties in the values of the energies can be neglected in all cases. 
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TABLE III: The contributions to the scalar polarizability for the 7pi/2 state in cesium. The corresponding 
energy differences and the absolute values of the lowest-order (DF) and final all-order electric-dipole reduced 
matrix elements are also listed. The energy differences are given in cm~^. Electric-dipole matrix elements are 
given in atomic units (eao), and polarizabilities are given in 10'' a^, where ao is Bohr radius. 



Contribution 




2;DF 


2-SD 


H '1^1/2 


Qo(7pi/2) 


_ main / ^ „N 
Q [TIS 1 


OS 


U.o / 


U.Zo 


-Zi ( DO 


n nnn 
-u.uuu 




{ S 


1 1 m 

11. Ul 


1 n 91 

lU.ol 


-OZOKJ 






8s 


9.53 


9.31 


2552 


2.487(35) 




9s 


2.04 


1.97 


5145 


0.055 




10s 


1.04 


1.00 


6535 


0.011 




lis 


0.68 


0.65 


7366 


0.004 




12s 


0.49 


0.48 


7904 


0.002 


„ tail c \ 

a [ns) 










0.012(12) 




5^3/2 


4.04 


1.52 


-7266 


-0.023 




6^3/2 


19.62 


17.99 


824 


28.74(70) 




7^3/2 


4.03 


6.56 


4283 


0.734(5) 




8^3/2 


2.39 


3.16 


6046 


0.121 




9^3/2 


1.63 


2.00 


7063 


0.042 




10d3/2 


1.21 


1.44 


7703 


0.020 


Q*-'('^d3/2) 










0.080(80) 


Total 










29.89(70) 



Next, we consider the examples of the polarizability calculation; one case is considered in detail for 
each of the nS, Npi/2, and Nd^i2 sequences of states. We consider the following sample 

cases: 9s, lpi/2i ^Vz|2^ ^^zjii 8c?5/2. In addition, we consider the 5(i3/2 and 6^5/2 calculations 
separately as they do not follow the pattern of all other Nd state calculations. These are also the only 
cases where the tail contribution is significant and represent interesting exception among the states 
that we have considered. 

We consider the 9s case first. The detailed breakdown of the 9s polarizability calculation is given 
in Table [TTl Each contribution to the main term, i.e. the contributions from the 6p, 7p, 8p, 9p, lOp, lip 
and 12p states are given separately, and the tail terms are grouped together for the npxji and np3/2 
contributions. The corresponding main term energy differences and the absolute values of the lowest- 
order (DF) and final all-order electric-dipole reduced matrix elements are also listed. The lowest-order 
values are given to illustrate the size of the correlation corrections for these transitions. The energy 
differences are given in cm~^. Electric-dipole matrix elements are given in atomic units (eao), f^nd 
polarizabilities are given in 10^ Oq. The core contribution is negligible in this case (0.015 in the units 
of Table nil and is not listed. We find that two of the transitions, 9s — 9pi/2 and 9s — 9p3/2, give 
dominant contributions while two other, 9s — 8pi/2 and 9s — 8^3/2, are large and have to be calculated 
accurately. We note that there is rather significant cancellation between the 9s — 9j>j and 9s — Spj 
contributions. The dominant contribution is this case may have been easily predicted simply based 
on the size of the energy intervals listed in the fifth column of the table. We also find that all other 
contributions with the exception of the 9s — lOpi/2 and 9s — IOP3/2 contribution are very small and 
may be simply omitted without the loss of accuracy. The main uncertainty comes from the uncertainty 
in the 9s — 9p3/2 transition. The precision our calculation in this case is expected to be very high as 
the correlation correction is small as illustrated by the comparison of the lowest-order and final values 
of the electric-dipole matrix elements. The final uncertainty is evaluated to be 0.7%. The breakdown 
of the calculation of the other N s polarizabilites considered in this work is similar to the one for the 
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TABLE IV: The contributions to the scalar and tensor polarizabilities for the IP3/2 state in cesium. The 
corresponding energy differences and the absolute values of the lowest-order (DF) and final all-order electric- 
dipole reduced matrix elements are also listed. The energy differences are given in cm~^. Electric-dipole matrix 
elements are given in atomic units (eao), and polarizabilities are given in 10"^ ajj, where ao is Bohr radius. 



Contribution 








E/i - Ejp^^^ 


ao(7p3/2) 


Q2(7p3/2) 


„ main /„ „\ 

a [ns] 


6s 


0.69 


0.59 


-21946 


-0.001 


0.001 




7s 


15.35 


14.32 


-3411 


-2.199(4) 


2.199(4) 




OS 


1 A 00 
14. Z5 


14. U 1 


Zol L 


6.\}0\6 ) 


-o.Uo(o ) 




9s 


2.86 


2.73 


4964 


0.055 


-0.055 




10s 


1.44 


1.38 


6354 


0.011 


-0.011 




Us 


0.93 


0.89 


7185 


0.004 


-0.004 




12s 


0.68 


0.65 


7722 


0.002 


-0.002 


„ tail r \ 

cx [ ns 1 










0.01(1) 


-0.01(1) 


a--"(nd3/2) 


54/2 


1.69 


0.58 


-7447 


-0.002 


-0.0013 




04/2 


O.OD 


O.U / 


D4z 


Q 71 / 1 0\ 
0.1 i(W) 


Z.9 1(0) 




74 /o 


2.11 


3.32 


4102 


0.098(1) 


0.079(1) 




84/2 


1.19 


1.54 


5865 


0.015 


0.012 




94/2 


0.79 


0.96 


6882 


0.005 


0.004 




104/2 


0.58 


0.68 


7522 


0.002 


0.002 


a*-'(nrf3/2) 










0.009(9) 


0.007(7) 




54/2 


5.02 


1.87 


-7350 


-0.017 


0.004 




64/2 


26.61 


24.35 


685 


31.6(7) 


-6.33(15) 




74/2 


6.30 


9.64 


4122 


0.825(6) 


-0.165(1) 




84/2 


3.55 


4.52 


5877 


0.127 


-0.025 




94/2 


2.37 


2.83 


6889 


0.042 


-0.009 




104/2 


1.75 


2.02 


7527 


0.020 


-0.004 












0.08(8) 


-0.02(2) 


Total 










37.52(75) 


-4.41(17) 



9s state with the exception of the 6s state. For all other cases, the dominant contributions come from 
the Ns — Npi/2 and the Ns — Np^/2 niatrix elements, while the other important contributions come 
from the Ns — {N — l)pi/2 and the Ns — {N — l)pz/2 matrix elements. The polarizability of the 6s 
state is overwhelmingly dominated by the contribution of the 6s — Qpi/2 and the 6s — 6^3/2 transitions. 
These two transitions add coherently and account for the 96% of the total value. The calculation of 
the 8s polarizability is described in detail in Ref. [3£||. We limited this work by the 12s state as the 
Hpj states needed for the calculation of the 13s polarizability do not quite fit inside of our cavity and 
the basis set energies of the 13p states deviate from the DF energies. 

The breakdown of the contributions to the lpi/2 and 7pz/2 polarizabilities is given in Tables Hill 
and IIV[ respectively. All tables illustrating the contributions to polarizabilities are structured in the 
same way. In the case of the 7pi/2 polarizability, the dominant contribution comes from a single 
transitions, 7pi/2 — 6^3/2, as none of the other levels are as close to the 7pi/2 levels as the 6^3/2 
level. The contribution from the next transition, 7pi/2 — 7^3/2, is significantly smaller, only 2% of the 
dominant contribution. Interestingly, the contributions of the 7pi/2 — 7s and 7pi/2 — 8s transitions, 
while being 10% of the main contribution, cancel out nearly exactly. We note that while significant 
cancellation is present for all other Npi/2 cases, it is the most severe in the case of the 7pi/2 state. 
The tail contribution is larger than for the Ns calculation but is still very small, 0.3%. We assume 
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TABLE V: The contributions to the scalar and tensor polarizabilities for the 8^3/2 state in cesium. The 
corresponding energy differences and the absolute values of the lowest-order (DF) and final all-order electric- 
dipole reduced matrix elements are also listed. The energy differences are given in cm~^. Electric-dipole matrix 
elements are given in atomic units (eao), and polarizabilities are given in 10'^ sq, where ao is Bohr radius. 



Contribution 


/3 


yDF 


ySD 


Ef3 - Esds^2 


00(86(3/2) 


02(86(3/2) 


(npi/2) 


6pi/2 


1.11 


1.30 


-16633 


0.00 


0.00 




7Pl/2 


2.39 


3.16 


-6046 


-0.06 


0.06 




oPl/2 


0.00 


y.io 


-ZWZ 


-1.47(2) 


1.47(2) 




9jDi /o 


50.96 


49.29 


-174 


-510(3) 


510(3) 




lOpi/2 


20.43 


14.02 


916 


7.85 


-7.85 






5.84 


4.50 


1592 


0.46 


-0.46 




120i /o 

f -L / 


3.07 


2.44 


2041 


0.11 


-0.11 


Lt \'f'Jr'l/2J 










0.2(2) 


-0.2(2) 


_main/„„ \ 


7P3/2 


1 in 
1.19 


l.o4 


-OODO 


m 
-0.01 


m 
-0.01 






2.97 


4.71 


-2020 


-0.40 


-0.32 




9P3/2 


23.02 


22.13 


-130 


-138.3(9) 


-110.6(7) 




10p3/2 


8.43 


5.53 


942 


1.18 


0.95 




llP3/2 


2.49 


1.83 


1610 


0.08 


0.06 




12p3/2 


1.32 


1.00 


2053 


0.02 


0.01 


a*""(nP3/2) 










0.04(4) 


0.03(3) 




4/5/2 


2.34 


2.49 


-3339 


-0.07 


0.01 




5/5/2 


19.18 


26.06 


-840 


-29.6(4) 


5.9(1) 




6/5/2 


70.91 


65.22 


518 


300(3) 


-60.0(6) 




7/5/2 


8.74 


0.33 


1337 


0.00 


0.00 




8/5/2 


5.71 


1.28 


1868 


0.03 


-0.01 


a*""(n/5/2) 










1(1) 


-0.2(2) 


Total 










-369(5) 


339(4) 



100% uncertainty in the tail contributions in all of our calculations for consistency. It is still negligible 
for all of the cases with the exception of the 5d calculation. 

As noted above, there are three types of the transitions contributing to the polarizabilities of the 
np3/2 states. The dominant contribution comes from the single transition as in the case of the 7pi/2 
polarizabilites, 7p^/2 — The contribution of the — 6(^3/2 transition is 10 times as small as 

the dominant one. Again, the contributions from the lp^/2 — 7s and 7p3/2 — 8s partially cancel, but 
the cancellation is not as complete as in the case of the 7pi/2 states. 

While the calculations of the scalar and tensor polarizabilities use the same matrix elements and 
energies and only differ by the angular factors, the uncertainty of the 7p3/2 tensor polarizability 
calculation (4%) is twice as high as that of the scalar polarizability owing to the significant cancellation 
of the terms contributing to the tensor polarizability. The relative accuracy of the calculation of the 
tensor polarizability calculation gradually improves to 1% for the IOP3/2 state but this uncertainty is 
still more than twice as high as the uncertainty of the corresponding scalar polarizability calculation 
(0.4%). The breakdown of all other npi/2 and np3/2 polarizabilites parallels the one of the 7pi/2 and 
7p3/2 states. 

The contributions to scalar and tensor polarizabilities for the 8(^3/2 and 8d^/2 states in cesium are 
given by Tables |V] and IVH For the 8^3/2 states, three contributions are dominant, 8^3/2 — 9j5i/2, 
8(^3/2 ~9pi/2, and 8c?3/2 —6/5/2 for both scalar and tensor polarizabilites. Unlike the case of the Np^/2 
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TABLE VI: The contributions to the scalar and tensor polarizabilities for the 8d5/2 state in cesium. The 
corresponding energy differences and the absolute values of the lowest-order (DF) and final all-order electric- 
dipole reduced matrix elements are also listed. The energy differences are given in cm~^. Electric-dipole matrix 
elements are given in atomic units (eao), and polarizabilities are given in 10"^ ag. 



(; Ti t V 1 r> 1 1 i" 1 ("1 Ti 


H 


yDF 

^0,Sd5/2 


ySD 


E/s Esd^/2 


rvn 1 r 1 n I 
HQ l^OU;5/2 ) 




a'"^-(nP3/2) 


6p3/2 


1.59 


1.81 


-16091 


0.00 


0.00 




7p3/2 


3.55 


4.52 


-5877 


-0.08 


0.08 




8p3/2 


8.82 


13.65 


-2031 


-2.24(2) 


2.24(2) 




9P3/2 


69.07 


66.57 


-141 


-765(5) 


765(5) 




IOP3/2 


25.43 


17.30 


931 


7.84 


-7.84 




IIP3/2 


7.51 


5.68 


1598 


0.49 


-0.49 




I2P3/2 


3.98 


3.10 


2041 


0.11 


-0.11 












0.2(2) 


-0.2(2) 


main / t~i \ 

a (^^^5/2 ) 


4.f5/2 


0.62 


0.67 


-3350.4 


0.00 


0.00 




5/5/2 


5.11 


6.85 


-851.3 


-1.34(2) 


-1.53(2) 




6/5/2 


18.97 


17.54 


506.6 


14.8(1) 


16.9(2) 




7/5/2 


2.37 


0.04 


1325.1 


0.00 


0.00 




8/5/2 


1.54 


0.42 


1855.9 


0.00 


0.00 


a'""(«^5/2) 










0.05(5) 


0.05(5) 




4/7/2 


2.79 


2.99 


-3350.7 


-0.06 


0.02 




5/7/2 


22.82 


30.60 


-851.6 


-26.8(4) 


9.6(1) 




6/7/2 


84.82 


78.43 


506.5 


296(3) 


-106(1) 




7/7/2 


10.61 


0.19 


1325.0 


0.00 


0.00 




8./7/2 


6.89 


1.86 


1855.8 


0.05(5) 


-0.02 












1(1) 


-0.4(4) 


Total 










-475(5) 


678(5) 



states, significant cancellations are observed between terms for both scalar and tensor polarizabilites. 
We would like to specifically note interesting problem with the 8(^3/2 — 7/5/2 transition. While the 
DF value for the transition is 8.74, the final all-order number is very small, 0.33 owing to extremely 
large correlation correction that essentially cancels the lowest order. We also note that the ah initio 
all-order value for this transition (0.73) significantly differs from the scaled values. While we assigned 
this value 100% uncertainty, the resulting uncertainty in the polarizability value is negligible. 

We observe similar problem with the 8(^5/2 — 7/7/2 transition as well as similar transitions for other 
values of N and n with the exception of the Sd — 4/ transitions. For the case of the 6(i — 5/ transition, 
the cancellation of the lowest order and the correlation correction is less severe. We note that the 
correlation correction to the previous transition in the sequence, 8(^3/2 — 6/5/2 is small, only 8%. 
Similar issue exists for the next in line transition, 8^3/2 — 8/5/2, but its contribution was too small to 
warrant its more accurate consideration. 

The 8(i3/2 — 7/5/2 and 8(^5/2 — 7/7/2 transitions are two of the very few transitions for which we 
conducted the scaling but did not list the values in the Table |T] of the recommended matrix elements 
as the uncertainties of these values are very high. In general, if the main term transition was not listed 
in Table U we used ah initio SD value and did not conduct the evaluation of the uncertainty. The 
contributions of these terms are small enough so their contribution to the total uncertainties would 
be negligible. Again, significant cancellations are observed between the terms. The polarizability 
calculation of the all other Nd is similar to the 8(i examples with the exception of the 5c? scalar 
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TABLE VII: The contributions to the scalar and tensor polarizabilities for the 5dg^2 and 6^5/2 states n cesium. 
The corresponding energy differences and the absolute values of the lowest-order (DF) and final all-order 
electric-dipole reduced matrix elements are also listed. The energy differences are given in cm~^. Electric- 
dipole matrix elements are given in atomic units (eao), and polarizabilities are given in 10^ Oq, where ao is 
Bohr radius. 



Contribution 


/3 


yDF 


ySD 






ao(5d3/2) 


a2(5d3/2) 


a'""(m/2) 


6pi/2 
7pi/2 


8.98 
4.04 


7.06 
1.52 




-3321 
7266 


-0.550(24) 
0.012 
0.002 


0.550(24) 
-0.012 
-0.002 


a'""(np3/2) 


6P3/2 
7p3/2 


4.06 
1.69 


3.19 
0.58 




-2767 
7447 


-0.134(6) 
0.002 
0.000 


-0.107(5) 
0.001 
0.000 


^■""'"(n/5/2) 


4/5/2 
5/5/2 
6/5/2 
7/5/2 
8/5/2 


10.66 
4.72 
2.90 
2.04 
1.55 


7.11 
3.34 
2.24 
1.66 
1.30 




9973 
12472 
13830 
14649 
15180 


0.186(27) 
0.033 
0.013 
0.007 
0.004 

0.059(59) 


0.037(5) 
-0.007 
-0.003 
-0.001 
-0.001 
-0.012(12) 


Total 












-0.352(69) 


0.370(28) 


Contribution 




yDF 

^/3,5d5/2 


ySD 

^/3,5d5/2 






00(5^5/2) 


02(5^5/2) 


«'""'"(np3/2) 

«*""(nP3/2) 


6P3/2 
7p3/2 


12.19 
5.02 


9.66 
1.87 




-2865 
7350 


-0.794(33) 
0.012 
0.002 


0.794(33) 
-0.012 
-0.002 


„ main / „ r \ 

a ('1/5/2) 
a'^"(n/5/2) 


4/5/2 
4/5/2 


2.84 
1.26 


1.93 
0.91 




9876 
12375 


0.009(3) 
0.002 
0.004(3) 


0.011(3) 
0.002 
0.004(3) 


a--"(n/7/2) 


4/7/2 
5/7/2 
6/7/2 
7/7/2 
8/7/2 


112.70 
5.64 
3.46 
2.44 
1.86 


8.62 
4.08 
2.73 
2.01 
1.57 




9875 
12375 
13733 
14551 
15082 


0.184(24) 
0.033 
0.013 
0.007 
0.004 

0.056(56) 


-0.066(9) 
-0.012 
-0.005 
-0.002 
-0.001 
-0.020(20) 


Total 












-0.453(70) 


0.691(40) 



polarizability calculation, which is anomalous and is discussed se par ately below. The calculation of 
the 7d, 9d, and lOd polarizabilities was discussed in detail in Ref. [23 |. 

The contributions to the scalar and tensor polarizabilities for the 5^3/2 and 5^5/2 states in cesium 
are given in Table IVIII We grouped small contributions of the 5d — npj and 5d — nf^ /2 transitions 
together with the tail in this table. Comparison of the 6^3/2 and 8^3/2 tables (as well as all the other 
nd^/2 contribution breakdowns) shows the 5^3/2 scalar polarizability case to be anomalous. In this 
case, none of the bd — nf energy denominators are small, and the largest contribution from n/5/2 states 
is still a third of the one from the dominant 6^3/2 — Qpi/2 transition. There is also no damping of the 



16 



TABLE VIII: Comparison of the Cs scalar polarizabilities with other theory and experiment. All values are 
given in 10^ Og. "Recommended value from Ref. [l^l, ''ab initio all-order value from Ref. [24l ]. °Ref. [3^. '^Re f.[^ , 
■^derived from the Ref. [T3 ] 7s — 6s Stark shift measurement and the 6s result from Q, ■'^Ref. [35| , ^Refs. (l^. |31|]. 
''derived from Ref. [ij] Dl line Stark shift measurement and the 6s result from [g, *Ref. [52l|. ■'derived from 
Ref. O D2 line Stark shift measurement and the 6s result from 8], ''Ref SSj, 'Ref. '541 "Ref. [§1 "Ref. 



State 


6s 


7s 


8s 


9s 


10s 


lis 


Present 


0.3984(7) 


6.238(41) 


38.27(28) 


153.7(1.0) 


478(3) 


1246(8) 


Ref. [31] 


0.394 


6.14 


37.9 


153 


475 


1240 


Theory 


0.3999(19)'' 


6.272'' 












0.4091^= 












Expt. 


0.4010(6)'' 


6.238(6)"= 


38.06(25)^ 




479(1)^^ 


1246(l)f 


State 


6pi/2 


7pi/2 


8pi/2 


9pi/2 


lOpi/2 


12s 


Present 


1.338(54) 


29.9(7) 


223(2) 


1021(7) 


3499(19) 


2866(30) 


Ref. [31] 


1.29 


29.4 


221 


1020 


3490 


2840 


Expt. 


1.3284(6)'' 


29.6(6)' 








2871(2)» 


State 


6p3/2 


7p3/2 


8p3/2 


9p3/2 


10p3/2 




Present 


1.648(56) 


37.5(8) 


284(3) 


1312(7) 


4522(19) 




Ref. [31J 


1.60 


36.9 


282 


1310 


4510 




Expt. 


1.641(2)^ 


37.9(8)'= 










State 


5^3/2 


6^3/2 


7^3/2 


8^3/2 


9d3/2 


10d3/2 


Present 


-0.352(69) 


-5.68(45) 


-66.7(1.7) 


-369(5) 


-1402(13) 


-4234(32) 


Ref. [31J 


-0.418 


-5.32 


-65.2 


-366 


-1400 


-4220 


Expt. 






-60(8)' 




-1450(120)'" 


-4185(4)" 


State 


54/2 


6^5/2 


7^5/2 


8^5/2 


9^5/2 


10d5/2 


Present 


-0.453(70) 


-8.37(55) 


-88.8(2.0) 


-475(5) 


-1777(14) 


-5316(38) 


Ref. [31J 


-0.518 


-7.95 


-87.1 


-472 


-1770 


-5300 


Expt. 






-76(8)' 




-2050(100)'" 


-5303(8)" 



remaining 5^3/2 ~ 1^/5/2 contributions observed for the higher 8d^/2 — ?^/5/2 transitions. Therefore, 
there is basis to assume that the DF tail is substantially overestimated. It may be overestimated by 
about 15-20% based on the comparison of the DF and the all-order matrix element values. As a result, 
the tail contribution is 25% of the total contribution of the 5^3/2 ~ sum and its uncertainty 

gives the dominant contribution to the uncertainty of the 5^3/2 scalar polarizability. We note that 
the 5(^3/2 — nf5/2 tensor polarizability tail is small with comparison to the dominant 5^3/2 ~ 
contribution, and its contribution to the total uncertainty is small. As a result, the 5^3/2 tensor 
polarizability calculation is similar to the 8^3/2 one. Its reduced accuracy is due to much larger 
correlation correction to the 8^3/2 — ()Pi/2 matrix element in comparison to the 8(^3/2 — one as 

illustrated by the comparison of the lowest-order and the all-order 5^3/2 — 6pi/2 and 8^3/2 — ^Pi/2 data. 
The analysis of the 5d^/2 polarizability is similar to that of the 5(^3/2 one. The main contribution to 
the uncertainty of the scalar polarizability comes from the 5ci5/2 — n-f7/2 tail and the uncertainties of 
the dominant terms are substantially larger than the uncertainties for the other Nd states for both 
scalar and tensor polarizabilites owing to large correlation correction of the corresponding transitions. 
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TABLE IX: Comparison of the Cs tensor polarizabilities with other theory and experiment. AU values are 
given in 10^ ag. "^Ref. ''Rei. [H, "Ref. [H "^Ref. [sj "Ref. ^ReL ^Ref. [H ''Ref. [H 



o uaue 


6p3/2 


7p3/2 


8p3/2 


9p3/2 


1 fin. 




Present 


-0.261(13) 


-4.41(17) 


-30.6(6) 


-135(2) 


-451(5) 




Ref. [31] 


-0.223 


-4.28 


-30.2 


-134 


-449 




Expt. 


-0.2624(15)" 


-4.43(12)** 


-30.7(1.2)"= 












-4.33(17)'* 














-4.00(8)"= 










State 


5^3/2 


6^3/2 


74/2 


84/2 


94/2 


104/2 


Present 


0.370(28) 


8.77(36) 


71.1(1.2) 


339(4) 


1189(10) 


3416(26) 


Ref. [31] 


0.380 


8.62 


70.4 


336 


1190 


3410 


Expt. 






74.5(2.0)^ 


333(16)"= 


1183(35)^ 


3401(4)3 


State 


5^5/2 


6^5/2 


74/2 


84/2 


94/2 


104/2 


Present 


0.691(40) 


17.33(50) 


142(2) 


678(5) 


2386(13) 


6869(34) 


Ref. [31] 


0.703 


17.00 


140 


675 


2380 


6850 


Expt. 






129(4)'' 


734(4)"= 


2660(140)"= 


6815(20)3 



7140(36)"= 



V. COMPARISON WITH OTHER THEORY AND EXPERIMENT 

Our final results for the scalar and tensor Cs polarizabilities and their uncertainties are compared 
with other theoretical and experimental values in Tables IVIIII and HXl respectively. As we noted above, 
the theory values for the 8s,7d,9d, and lOd polarizabilites from Refs. [12,[3^ differ very shghtly from 
the present values since they are obtained using the same values of the matrix elements but more 
accurate energies. Therefore, we do not quote theory values from Refs. [l^iHEl separately in Tables lVIIII 
and IIXI The experimental values for the 7s, 6p^ /2, and 6p^/2 states are obtained by combining the 
most accurate measurements of the 7s — 6s 6pi/2 — 6s jllj . and 6^3/2 — 6s [l^ Stark shifts with 
the recent measurement of the 6s polarizability [S'l, respectively. We find excellent agreement of our 
values with high-precision measurements of Refs. ^8,, ,10, 11, 12, 21, 35, 56] . Disagreements with older 
values for the Nd states are discussed in detail in Ref. [2^ . In all cases where the new measurements 
are available, our data support most precise measurements. In particular, we find that our method 
works very well for even such highly-excited states as 12s and IQd. 

We also compare our values with the van Wijngaarden and Li [sij work where the extensive calcula- 
tions of the polarizabilities of cesium Ns, Np, Nd, and Nf states were carried out using the Coulomb 
approximation. Our values are in excellent agreement with those results for higher excited states where 
the method of Ref. 'sij is expected to work well. 

VI. CONCLUSION 

We have carried out a systematic study of the Cs electric-dipole static polarizabilities for the Ns 
{N = 6 — 12), Npj (N = 6 — 10), and Ndj {N = 5 — 10) states using the relativistic all-order method. 
The recommended values for the polarizabilites of all these states are given and their uncertainties are 
estimated. This work involved the calculation of 317 electric-dipole transition in Cs. Recommended 
values for the 91 transitions that give the dominant contributions to the polarizabilities are presented 
together with their uncertainties. Our polarizability values are compared with other theory and exper- 
iment. Our data are found to be in excellent agreement with the high-precision measurements. These 
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calculations provide a theoretical benchmark for a large number of Cs clcctric-dipole matrix elements 
and polarizabilities. 
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